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Quantum error correction (QEC) is an essential element of physical quantum information pro- 
cessing systems. Most QEC efforts focus on extending classical error correction schemes to the 
quantum regime. The input to a noisy system is embedded in a coded subspace, and error recovery 
is performed via an operation designed to perfectly correct for a set of errors, presumably a large 
subset of the physical noise process. In this paper, we examine the choice of recovery operation. 
Rather than seeking perfect correction on a subset of errors, we seek a recovery operation to max- 
imize the entanglement fidelity for a given input state and noise model. In this way, the recovery 
operation is optimum for the given encoding and noise process. This optimization is shown to be 
calculable via a semideflnite program (SDP), a well-established form of convex optimization with 
efficient algorithms for its solution. The error recovery operation may also be interpreted as a com- 
bining operation following a quantum spreading channel, thus providing a quantum analogy to the 
classical diversity combining operation. 

PACS numbers: 03.67.Pp, 02.60.Pn 
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I. INTRODUCTION 

Any implementation of quantum computing or com- 
munications requires a strategy for error mitigation. 
Indeed, the development of quantum error correction 
(QEC) schemes was an important early step in moving 
quantum computing from an interesting theoretical idea 
to an exciting field with potential for ground-breaking 
technological implementations 1]. The importance of ef- 
ficient and optimum error mitigation only increases as 
the field advances. 

The earliest efforts in QEC used encoding techniques 
modified from classical error correction schemes 1, 13, 0, 
0, . Further analysis 0, laid the foundation for QEC 
theory, noting that the important metric is how faithful 
the statistics of the corrected state remain to the ideal 
behavior. While that observation suggests quantum error 
mitigation is thus an optimization problem, most of the 
subsequent work in the field has appropriately focused on 
perfect recovery from a set of errors. This emphasis has 
allowed many techniques to be borrowed from classical 
error correction and enabled important feasibility studies 
in quantum computing. It is not, however, the only way 
to consider controlling for quantum errors 8j. 

Recently, some authors 0, ^^^^ returned to exam- 
ining quantum error mitigation as an optimization prob- 
lem. The essential properties of a quantum state are the 
statistics of any observable outcome; these are completely 
encapsulated in the density operator p of the state. Noise 
is introduced by the operation £ which can be thought of 
as a noisy quantum communications channel. Thus, the 
goal of any error correction scheme is to design a recovery 
operation IZ such that the recovered state is as faithful 
a representation as possible of the input, judged by how 
well the statistics of observables are preserved. The opti- 



mum recovery minimizes the 'distance' between an input 
density p and the output 1Z{£{p)). This operation may 
differ from the more traditional QEC recovery operation; 
such differences illustrate further the contrast between 
quantum and classical error correction. To distinguish 
this approach from QEC, we use the term quantum error 
recovery (QER). It should be emphasized that the opti- 
mum QER recovery operation is dependent on a given 
input density, encoding operation, and noise model. 

The paper is organized as follows. In Sec. ^ we define 
the parameters for optimum QER. Section imi derives a 
representation of a quantum channel based on a single 
positive operator. In section HVl the optimum recovery 
operation is cast as a semidefinite program. Section 
interprets the recovery operation as an optimum combin- 
ing problem, and illustrates the computational benefit of 
such an interpretation. In Sec. IVIl QER operations are 
derived for the amplitude damping channel using codes 
encoding one qubit into four and five qubits. 



II. OPTIMUM QER 

Most QEC procedures are designed on the principle of 
'perfect' correction of arbitrary single qubit errors. Such 
a design postulates that single qubit errors are the dom- 
inant terms in the noise process; thus a scheme that cor- 
rects arbitrary single qubit errors and ignores higher or- 
der terms will sufficiently mitigate the noise. Pursuit of 
this approach has led to important results on the feasi- 
bility of quantum error correction, and indeed is a rea- 
sonable model for noise processes accurately described by 
the lower order terms. However, one may reasonably ask 
how well this generic approach succeeds in specific cases. 

With every quantum code, in the current paradigm. 
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there is an associated recovery operation designed to per- 
fectly correct the dominant errors. This 'traditional' re- 
covery (referred hereafter as the QEC recovery operation) 
applies a syndrome measurement to determine which er- 
ror occurred, and a correction operation dependent on 
the observed syndrome. For a given code and error pro- 
cess, the QEC recovery operation may not provide the 
most effective safeguard from error. Depending upon 
the form of the error process, an alternate recovery op- 
eration may be designed that better preserves the in- 
put state, based on some measure of statistical 'close- 
ness' between the input density p and the output density 
lZ{£{p)). Commonly used metrics^ for quantum infor- 
mation arise from the fidelity^ defined to be 

a) = tr VpV2^pi/2 (1) 

where p and a are density operators. The fidelity F takes 
a value between and 1, where 1 indicates that the states 
are identical. Using Fq to represent any fidelity-based 
measure of channel performance, the error recovery opti- 
mization problem is to find 

IV = arg m^x (7^ o f ) , (2) 

where the maximization is over all valid quantum opera- 
tions. 

We vaguely described the measure of channel perfor- 
mance, Fo, as fidelity-based. In fact, the specific choice 
of Fo infiuences the feasibility of the optimization, as 
well as the interpretation of the result. There are three 
main choices of fidelity-based measures of performance, 
the minimum fidelity, the ensemble average fidelity, and 
the entanglement fidelity 11]. 

The minimum fidelity has the advantage of bounding 
the performance by seeking the worst case input state. 
The metric Fq is the minimization over all inputs of the 
input-output fidelity. The optimization becomes two- 
fold:2 

1V = argmaxminF(p,7^(f (p))), (3) 

71 p 

where F is the fidelity given in By virtue of the 

minimization over p, one need not assume anything about 
the input state. This was the metric of choice in 
first establishing the theory of QEC. The disadvantage 
arises through the complexity of the metric; indeed the 
optimization problem is now over two sets, and the metric 
Fo is not linear with respect to 7Z. Efficient routines that 
have been developed for solving the optimization problem 
of ^ are sub-optimum |lOj . 



We obtain a more tractable optimization problem if we 
are willing to assume some form for the input distribution 
p, particularly if the metric Fq can be reasonably defined 
as linear in the recovery operation TZ. While assuming a 
form for p makes the solution less general, it illuminates 
important characteristics of quantum error mitigation by 
enabling construction of an optimum recovery operation. 
Given a value of p, we may use either the ensemble aver- 
age fidelity or the entanglement fidelity. 

Ensemble average fidelity models the input as being 
in a state pi with probability Pi. We define an arbi- 
trary quantum channel B on the Hilbert space H as 
B : jC{H) ^ where C{H) indicates the set of 

bounded linear operators on H. The measure by which B 
preserves the input is then the average squared fidelity: 

F{pi,Pi,B)=Y,P^PiP^,BiP^)?■ (4) 

F is linear in B when each pi describes a pure state; for 
linearity we must assume more than just the density of 
the input. 

Entanglement fidelity 12] arises from the mathematical 
concept of mixed state purification. Any mixed quantum 
state can be represented as a subsystem of a pure state 
in a larger Hilbert space. The subsystem is mixed due to 
the entangled nature of the pure state. Consider a mixed 
state p e OH). By defining a reference space A, we may 
denote p as a subsystem of a pure state: 

p = iXA \AH) {AH\ (5) 

where \AH) is a pure state in the space A^7i. The entan- 
glement fidelity then measures how faithfully B maintains 
the purification (or equivalently, how well it preserves the 
entanglement). It is given by^ 

F{p, B) = {AH\ Ia ® B{\AH) {AH\) \AH) , (6) 

the squared fidelity of the input \AH) with the output of 
the channel 2"^ B. 

With entanglement fidelity as Fo in Q, we define the 
optimum recovery operation for the error process S and 
input distribution p as 

= argmaxF(p,7^of). (7) 



III. CPTP MAPS AND POSITIVE OPERATORS 

A valid quantum operation must be completely pos- 
itive (CP) and trace-preserving (TP) This require- 
ment follows directly from the postulates of quantum me- 
chanics wherein the evolution of a closed quantum system 



^ While not technically a metric, the fidelity is a useful measure of 

performance for quantum states. 
^ For simplicity of notation, the minimization is shown over all 

densities p. In fact, it is sufficient to minimize over pure state 

inputs IJJ . 



^ We use F to signify both the fidelity and the entanglement fi- 
delity. The distinction should be obvious from context, as the 
arguments for the fidelity are two density operators, whereas for 
the entanglement fidelity, the arguments are a density operator 
and a channel mapping. 



3 



is unitary. Let 5 be a CPTP map from C{H) ^ £(/C). 
The most familiar representation of the CPTP map is 
the Kraus (or operator sum) form, where the mapping is 
specified by a set of operators {E^} known as the oper- 
ator elements jl3ij. The channel output is given by 

Sip) = ^EkpEl (8) 

k 

and the CPTP constraint is met when 

J2ElEk=I, (9) 

k 

where / is the identity operation on C{TL). 

While a properly constrained set of operators fully 
specify a quantum channel, the Kraus form is an incon- 
venient one for optimization. The most obvious inconve- 
nience is the many-to-one correspondence between sets of 
operator elements and channel mappings. For this rea- 
son, we will utilize the representation of a channel map- 
ping by a positive semidefinite operator on £(/C H), 
which correspondence is one-to-one [3 • We will refer to 
the operator description of CPTP maps as the superop- 
erator form 15, 16]. 

The superoperator form may be derived by recogniz- 
ing that the space of bounded linear operators forms a 
Hilbert space. It is convenient to have a general method 
to convert the operator to a ket notation. Let C be a 
bounded linear operator from 7^2 to TYi: C e C{TL2^Hi). 
We define a ket in the Hilbert space TYi associated 
with C as^ 

\C)) = Y.''Mj)2. (10) 

where and {|j)2} are orthonormal bases for Hi 

and 7^2, respectively, and Cij = i{i\C\j)2 is the matrix 
element of C on these bases. 

Two useful relations follow directly from this defini- 
tion. The first one, 

A^B\C)) = \ACB^)), (11) 

applies whenever the dimensions of A, and C indicate 
that ACB^ is a valid operator.^ The second relation 
applies for Ci,C2 G C{H2^Hi): 

tTnA\Ci)){{C2\] = C,CleC{n,). (12) 

From these relations, we see that the channel mapping 



^ The notation |-)} is used to emphasize that these kets represents 

operators, not quantum states. 
^ Note the symbol ^ indicates the transpose with respect to the 

specified basis without conjugation. 



E : jC{H) ^ >C(/C) can be given by 

£{p) = Y.^'^p^i 

k 

= J2^rn[\Ekp)){{E,\] 

k 

= ^tr„ [l®p^\Ek)){{Ek\] 

k 

= trn[l®P^Xs], (13) 

where Xg = \Ek)){{Ek\. The trace-preserving prop- 
erty that 

tijcSip) = l = tTn [p^tijc [Xs]] (14) 

for all density operators p G JC{H) can be stated as 

tTjcXs = I eCin). (15) 

In the superoperator form, the entire mapping £ is spec- 
ified by the positive operator Xg. 

IV. OPTIMUM RECOVERY VIA 
SEMIDEFINITE PROGRAMMING 

The problem given by Q is a convex optimization 
problem and we may approach it with sophisticated tools. 
Particularly powerful is the semidefinite program (SDP) 
17], where the objective function is linear in an input 
constrained to a semidefinite cone. Indeed, the power of 
the SDP was a primary motivation in choosing to max- 
imize the entanglement fidelity, which is linear in the 
quantum operation 1Z. 

The definition of entanglement fidelity given in Q is 
intuitively useful, but awkward for calculations. An eas- 
ier form arises when operator elements {Bi} for B are 
given. The entanglement fidelity is then 

F{p,B) = J2MpBi)f. (16) 

From ([TB|) . we may derive a calculation rule for the 
entanglement fidelity when the channel B is expressed 
as in the superoperator form. Using ([T^ . we see that 
tr^^p = tT\Bi)){{p\ = {{p\Bi)). Inserting this into (|TB|) . we 
obtain the entanglement fidelity in terms of Xs- 

F{p,B) = J2{{pm){{BM 

i 

= {{p\Xb\p)). (17) 

Armed with this expression for the entanglement fi- 
delity, we may now express {TJ in a form readily seen to 
be a semidefinite program. To do this, we must consider 
the form of the composite channel 1Zo £ : OH) ^ >C(7Y) 
expressed as a positive operator on 7i ^ 7i. If the oper- 
ator elements for each channel are {Ri} and {Ej}, then 
the operator X-jzs is given by 

ij 



Applying (|TT|) . this becomes 



€ ^ c c{K) 



n{8{pi'')) €CC C{K) 



Xtzs = ^/®£^f|i?i))((i?i|/®£;* 

= ® Ef)Xn{I(E) E*), (19) 

where the * represents complex conjugation, without 
transposition. The entanglement fidelity is then 

F{p,no£) = Y,{{p\{I®Ej)Xn{I®E*)\p)) 

3 

= trXnCp^s (20) 

where 

Cp,e = Y.I®E*\p)){{p\I®Ej 

3 

= Y,\pE])){{pE]\. (21) 

3 

We may now express the optimization problem Q in 
the simple form 

■^Up — argniaxtr(XCp,£:) 

such that X > 0, tr-^X = /. (22) 

This form illustrates plainly the linearity of the objective 
function and the semidefinite and equality structure of 
the constraints. Indeed, this is the exact form of the opti- 
mization problem in ^3 which first pointed out the value 
of semidefinite programs (SDP) for optimizing quantum 
channels. 

The value of an SDP for optimization is two-fold. 
First, an SDP is a sub-class of convex optimization, and 
thus a local optimum is guaranteed to be a global op- 
timum. Second, there are efficient and well-understood 
algorithms for computing the optimum of a semidefinite 
program. These algorithms are sufficiently mature to be 
widely available. Thus, by expressing the optimum re- 
covery channel as an SDP, we have explicit means to 
compute the solution numerically. 

V. QUANTUM DIVERSITY COMBINING 

In the preceding analysis, the method of encoding has 
been implied by the choice of p. Indeed, in most treat- 
ments of QEC the input density is restricted to a sub- 
space called the quantum error correcting code (QECC). 
If Pc is a projection operator onto the code subspace, 
then p = PcpPc implies that the state p is within the 
code subspace. The channel is typically defined such that 
the input and output spaces are identical {i.e. H = JC) 
and the noise process generally perturbs the system from 
the code subspace. While this representation is a per- 
fectly legitimate model for the error process, and conve- 
nient when viewing QEC in a mode comparable to clas- 
sical error correction, the dimensionality is unnecessarily 
high for the optimization routine. 
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FIG. 1: Transform from standard quantum error correction to 
quantum spreading channel. By considering the error channel 
8 together with the encoding Uc and the recovery operation 
IZ with the decoder t/^, the dimension of the SDP may be 
reduced by a factor of (dim/C/ dimTY)^. Note that the output 
of the encoder Uc and the recovery channel IZ live on the code 
subspace C. 



Recall that X^z is a Hermitian operator on the space 
JC H. The optimization thus has d^^d^ real parame- 
ters. Even for a [[5^ ^^.^1]^ code, the smallest for arbitrary 
single qubit errors [71 [l9^. this optimization then has 2^^ 
optimization variables. The high dimensionality can be 
alleviated somewhat by embedding the encoding into the 
noise process, and redefining ^ as a quantum spreading 
channel where dn < d)c- 

The transform to a quantum spreading channel is illus- 
trated in Fig. nj Consider the noise process 8 : JC JC 
with operator elements {Ei}^ input density p G JC{JC), 
and code projector 

Pc = J2 \^)LL{n\ (23) 

n 

where |n)^ G JC are the logical states of the code. Since 
the input is in the codespace, p is preserved by the code 
projector: p = PcpPc- We can reduce the dimensional- 
ity of the optimization by transforming this problem. Let 
H he di dn dimensional Hilbert space with orthonormal 
basis {\n)^}. H can be considered the space of the in- 
formation source. The encoding process is accomplished 
by the operator Uc = Ylt^ that maps the basis 

states of Ti to the logical states in /C. Uc is an isometry: 
uIjUc = P Note that 

UcUl = ^|n)ii(n| 

n 

= Pc (24) 



The notation [[n,/c,(i]] refers to a quantum code encoding k 
qubits of information into a n qubit system. The parameter 
d refers to the weight of the code , 1 Li . 
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If we redefine p' = U^pUc ^ ^{^) and the operator 
elements E[ = EiUc^ then we see the processes 8{p) and 
£'{p') are identical: 

£'{p') = Y.EiUcUl;pUcUl;E\ 
= Y^EiPcpPcEj 

i 

= £{P)- (25) 

By enacting such a transformation, the optimization di- 
mension is reduced from d!^ to d^d^. For the [[5,1,3]] 
code, the reduction is from 2^^ to 2^^. 

The above transformation illustrates an alternative in- 
terpretation of recovering an encoded quantum state af- 
ter an error process. We may instead consider an unen- 
coded state input into a quantum spreading channel^ a 
channel in which the output dimension is greater than 
the input dimension {i.e. dim/C > dimTY). The recovery 
operation is an attempt to combine the spread output 
back into the input space, presumably with the intent 
to minimize information loss. The recovery operation is 
then the quantum analog to the classical communications 
concept of diversity combining. 

Classical diversity combining describes a broad class 
of problems in communications and radar systems. In its 
most general form, we may consider any class of trans- 
mission problems for which the receiver observes multiple 
transmission channels. These channels could arise due 
to multi-path scattering, frequency diversity (high band- 
width transmissions where channel response varies with 
frequency), spatial diversity (free-space propagation to 
multiple physically separated antennas), time diversity, 
or some combination of the four. Diversity combining is a 
catch-all phrase for the process of exploiting the multiple 
channel outputs to improve the quality of transmission 
{i.e. by reducing error or increasing data throughput). 

On its face, an extension of diversity combining to the 
quantum communications regime is unclear. Indeed, the 
simplest way to understand classical diversity - receiv- 
ing multiple copies of the input, independently affected 
by the channel - might also lead one to conclude that 
the concept is not applicable to quantum communica- 
tion. After all, the fundamental postulates of quantum 
mechanics lead to the "no-cloning" theorem; multiple, 
separable copies of the input state are ruled out by the- 
orem. A more careful consideration of classical diversity 
combining, however, illuminates the parallels between 
the quantum and classical viewpoints. 

In a general description of classical diversity, the input 
signal is coupled through the channel to a receiver sys- 
tem of higher dimensionality. Consider a communication 
signal with a single input antenna and N receive anten- 
nae. Often, the desired output is a signal of the same 
dimension as the input, a scalar in this case. Diversity 
combining is then the process of extracting the maximum 
information from the A/'-dimensional received system. In 
most communications systems, this combining is done at 



either the analog (leading to beam-forming or multi-user 
detection) or digital (making the diversity system a kind 
of repeater code) levels. Thus, the natural inclination is 
to equate diversity combining with either beam-forming 
or repeater codes. The most general picture, however, is 
that of recombining the received signal that was spread 
by the channel from the input system to a higher dimen- 
sional output system. Thus, it is appropriate to consider 
a quantum spreading channel to be quantum diversity 
channel, and the recovery operation is a quantum diver- 
sity combiner. 



VI. EXAMPLES 

To illustrate the potential benefit of optimum QER, 
we numerically compare the performance of the opti- 
mum QER with QEC recovery for two encoding schemes. 
First, we examine the [[5, 1, 3]] code in the presence of the 
amplitude damping channel. Second, we compare the 4- 
bit amplitude damping code from 8] with the optimum 
QER. The latter is a particularly apt choice for compar- 
ison, as the code was designed for a specific channel and 
sought to only approximately satisfy the quantum error 
correcting conditions. 

In both cases, the noise channel £a is the amplitude 
damping channel, which for a single qubit has operator 
elements 



^0 



1 

yr^ 



7 



and El 









(26) 



This channel is a commonly encountered model, where 
the parameter 7 indicates the probability of decaying 
from state |1) to |0) {i.e. the probability of losing a 
photon). Amplitude damping is a logical choice to il- 
lustrate the benefits of optimum QER as the operation 
is not symmetric with respect to |0) and 



A. Five-Qubit Code 

The five-qubit code was independently discovered by 
1^ and 0. We will here follow the treatment in 11], and 
specify the code via the generators {^1, ^2, ^3, ^4} and the 
logical Z and X operations specified in Table IVOl The 
code subspace C is the two-dimensional subspace that is 
the +1 eigenspace of the generators gi. The logical states 
|0)l and are the +1 and —1 eigenkets of Z on C. 

To compute the optimum recovery for this code, we 
assume that the logical states are equally likely, that 
is p = ||0)ll(0| + ||1)ll(1|. (This choice of p in fact 
assumes nothing about the choice of codewords; rather 
it is the maximum entropy distribution constrained to 
the code space C.) The noise channel is f the am- 
plitude damping channel acting independently on each 
qubit. For each choice of the parameter 7, the optimum 
recovery operation 7^^ is computing according to Q. We 
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Name 


Operator 


91 


XZZXI 


92 


IXZZX 


93 


XIXZZ 


94 


ZXIXZ 


Z 


zzzzz 


X 


xxxxx 



TABLE I: Generators and logical operations of the five qubit 
{i.e. [[5, 1, 3]]) code. 



compare the entanglement fidelity F{p^7lo6f^) for both 
IZ^p and T^QEC in fig. El 

Figure El illustrates clearly the difference between op- 
timum QER and QEC recoveries for large values of 7. 
It is also instructive to compare the techniques for small 
values of 7. We do this numerically by calculating the 
polynomial expansion of IZ o £) as 7 goes to zero. 
The entanglement fidelity for the optimum QER has the 
form no£)^l- 1.1667^ + 0(7^)- In contrast, the 
QEC recovery is noS)^:^!- 2.67^ + 0(7^). 




code. The code only approximately meets the QEC con- 
ditions, and as a result, the 'corrected' state is somewhat 
distorted from the input, even when the dominant error 
term occurs. In this way, the procedure is based upon 
principles similar to the optimum QER we have devel- 
oped. 

The logical states are given by 



1 



|0)l = ^(10000) + 11111)) 
|1)l = ^(10011) + 11100)), 



(27) 
(28) 



and the recovery operation is specified by the circuits in 
Figure 2 of 8] . We note that the recovery operation de- 
pends explicitly on the parameter 7. We compare the re- 
covery of Leung et. al. with the optimum QER computed 
according to {TJ, once again assuming the completely 
mixed input density p = ^|0)ll(0| + ^|1)ll(1|- The nu- 
merical comparison, for various values of 7 is provided in 
fig. IHl As 7 goes to zero, the entanglement fidelity for 
the optimum QER is numerically determined to have the 
form F(p, lloS) ^l- 1.257^ + 0(7^). In contrast, the 
Leung et al recovery is F(p, T^of) ^ 1-2.757^ + 0(7^). 
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FIG. 2: Entanglement fidelity vs. 7 for the 5 qubit code and 
the amplitude damping channel Sf^. 7 refers to the damp- 
ing parameter of the channel. Entanglement fidelity with no 
recovery operation is included for comparative purposes. 



FIG. 3: Entanglement fidelity vs. 7 for the 4 qubit amplitude 
damping code of Leung et. al. ^ . 



B. Four-Qubit 'Approximate' Code 

In 0, Leung et. al. recognized the advantage that 
channel- specific error recovery schemes can have over 
generic QEC routines. To illustrate the advantage, they 
designed a code specific to the amplitude damping chan- 
nel with good performance for small 7 that only required 
4 qubits of overhead, in contrast to the generic five qubit 



C. Commentary on Examples 

It is not surprising that the optimum QER opera- 
tion outperforms the 5 qubit code or the Leung code. 
It is, however, noteworthy that the difference can be 
significant, for even relatively small values of 7. We 
may conclude that channel specific recoveries can signifi- 
cantly improve the performance of error correcting proce- 
dures. This conclusion was shared by Leung et. al, who 
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lamented the lack of a general method to design such 
channel-specific schemes. The SDP formalism outlined 
in this paper, provides such a general method. 

Perhaps the most notable improvement in optimum 
QER can be noted in the 'no recovery' comparisons in 
figs. El and [HI These curves represent the entanglement 
fidelity of a single qubit transmitted through the noisy 
channel. For values of 7 where the recovered entangle- 
ment fidelity lies below the 'no recovery,' we see that the 
error mitigation procedure does more harm than good. 
Performing optimum QER as opposed to QEC signifi- 
cantly extends the values of 7 for which error mitigation 
is valid. This suggests optimum QER will be a particu- 
larly valuable technique for noisier systems. 

Finally, it is worth noting the duality between opti- 
mum QER and optimum encoding. We have derived 
the optimum operation IZ for a given encoding and noise 
process. By the same process, one may derive the opti- 
mum encoding given a recovery operation and noise pro- 
cess. This can most easily be seen by noting that the 
encoding operation is a valid quantum operation, and in 
fact, a spreading operation; it is thus subject to the same 
semidefinite cone constraints as the recovery operation. 
In a manner similar to those suggested by |2ii] and 
one can conceivably obtain a channel-specific error re- 
covery scheme by alternatively holding the recovery fixed 
and optimizing the encoding, and holding the encoding 
fixed and optimizing the recovery. Full analysis of such 
a technique is deferred for future consideration. 



VII. CONCLUSION 

The structure of quantum operations allow quantum 
error correction to be approached as an optimization 
problem. Specifically, optimum recovery of an encoded 
quantum state from an error process can be solved nu- 
merically using semidefinite programming when optimal- 
ity is interpreted as a maximization of the entangle- 
ment fidelity. This analysis suggests the ability to sys- 
tematically search for recovery operations for compli- 
cated error schemes beyond those readily analyzed and 
corrected through more traditional QEC methodologies. 
This problem is in general, the optimum combining op- 
eration following a quantum spreading channel, and thus 
a quantum parallel to the classical problem of diversity 
combining. 
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